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Abstract 

We investigate the local spin polarization texture of Landau levels under Rashba spin-orbit 
coupling in bulk two-dimensional electron gas (2DEG) systems. In order to analyze the spin 
polarization as a function of two-dimensional coordinates within the 2DEG, we first solve the 
system eigenstates in the symmetric gauge. Our exact analytical wavefunction solutions are shown 
to be gauge invariant with solutions obtained in the commonly used Landau gauge. We illustrate 
the two-dimensional spatial spin profile for a single Landau level and suggest means to measure 
and utilize the local polarization in practice. 

PACS numbers: 72.25.Dc,71.70.Di,71.70.Ej 



1 



I. INTRODUCTION 



Spin dependent transport phenomena in low dimensional systems have attracted consider- 
able attention in recent years because of their potential application in information processing 
and storage devices [l-3|. A paradigmatic proposal is the spin field-effect-transistor which 
utilizes the gate-controllable [4] Rashba spin-orbit coupling (SOC) js, gI in two-dimensional 
electron gases (2DEGs) to control the spin rotation of electrons as they propagate across 
the device [t-O]. The Rashba SOC results from the structural inversion asymmetry of the 
microscopic confinement potential formed at the interface of semiconductor heterostructures 
5|, l6|. There is much interest recently in 2DEG systems with SOC and external magnetic 
fields. By applying a perpendicular magnetic field to the 2DEG system, the SOC competes 
with Zeeman spin-splitting and this interplay leads to further modification of band structure 
and other interesting results. A few examples are resonant spin-Hall conduct ance_^ 10|, 111] 



12| 



mod- 



due to induced degeneracies of Landau levels at certain values of magnetic field 
ified magneto-optical transition spectrums [3], beating patterns in the density of states 



15| which differs from the 



and longitudinal resistivity [l^, and altered Hall conductance 
quantized values in the integer quantum Hall (IQH) regime. We note that previous works 
(including all of the above) perform their analyses in the Landau gauge. While the use of the 
Landau gauge is perfectly valid due to gauge invariance, the form of the wavefunction in this 
gauge does not capture the natural, rotational symmetry of the eigenstates. For example, in 
the presence of Rashba SOC, the spin polarization of eigenstates exhibits interesting spatial 
textures whose features cannot adequately be reflected by the wavefunctions obtained in the 
Landau gauge. The study of the locally varying spin polarization within a 2DEG may have a 
number of interesting applications. For instance, a well-controlled spin texture with distinct 
spatial modulation may be used as a resolution test for surface spin probe techniques. Ad- 
ditionally, it may be possible, by means of some localized probes, to harness an efficient spin 
current source from spatial regions with high spin polarization. Here, the spatial separation 
of spins is reminiscent of the optical dispersion (spatial separation of optical frequencies) 
found in monochromators, suggesting that it may be used as a form of spin filter. 
In this article, we theoretically study the local spin polarization of Landau levels in the pres- 
ence of Rashba SOC within an infinite 2DEG. To do so, we first present analytical solutions 
of the eigenstates of the system in the rotationally isotropic symmetric gauge. We demon- 
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strate the gauge equivalence of our solutions with previously known solutions obtained in the 
Landau gauge. Finally, we show the spatial distribution of the spin components of Landau 
level states in the presence of Zeeman and Rashba SOC effects. 



II. THEORY 

A. Landau levels in the symmetric gauge 

We first solve the Landau level wavefunctions without the Rashba SOC and Zeeman inter- 
actions in the symmetric gauge i.e. the IQH states. The wavefunctions form an orthonormal 
set which we use as the basis functions in solving the complete system. Under an external 
vertical magnetic field B, the Hamiltonian of a spinless and otherwise free electron in a 
2DEG is written as I-Lq = n^/2me, where 11 = p + eA is the covariant momentum under 
the vector potential A which satisfies curl A = B, — e is the electron charge, and mg the 
effective electron mass. Fixing the magnetic field B does not uniquely determine the vector 
potential, i.e. there is a gauge freedom. For a magnetic field that is perpendicular to the 
plane of the 2DEG, pointing in the 5-direction by convention, the Landau (L) gauge is given 
by A'^ = {—BzU, 0, 0) whilst the symmetric (sym) gauge is given by A*^^™ = Bz/2{—y, +x, 0), 
where Bz is the magnetic flux density (in Tesla) of the external field and x, y are spatial co- 
ordinates in the plane of the 2DEG. In the presence of a uniform magnetic field, the system 



exhibits both translational and rotational symmetry about the ;z-axis. Under t 



le Landau 



20| 



gauge, it is well known that the solutions to the Hamiltonian are of the form 

^ni^^y) = exp{ik^x)-^n[{y -yo)/r)], (1) 



where j/o = hk^/eB^ is the y coordinate of the cyclotron center, r = ^Jh/eBz is the mag- 
netic length and ipn {n, an integer) are the normalized nth order Hermite polynomials. The 
wavefunction characterizes the nth discrete Landau level in the presence of a mag- 

netic field, with corresponding quantized energy spectrum En = f\uj{n + i). Although the 
choice of Landau gauge preserves the translational symmetry of the system, the rotational 
invariance is lost in Eq. ([T]). In describing the circular Landau orbits of electrons which 
form in the presence of B fields, it is more natural to use the rotationally isotropic sym- 
metric gauge. The use of the symmetric gauge has been applied previously to analyze other 
systems exhibiting rotational symmetry, e.g. in 2D two-electron systems [16j], and quantum 



see, for example, H, Q), and QDs in 



dots (QDs) in 2D parabolic confinement potentials 
radially symmetric hard- wall potentials with SOC 19|, in the presence of magnetic fields. 
Under this choice of gauge, it is convenient to define the complex variable z = x + iy to 



represent spatial coordinates within the 2DEG, and introduce the operators 



20|: 



In analogy to the harmonic oscillator, the Hamiltonian can be rewritten in terms of a and 
a\ i.e., T-Lo = fijjj{a)a + i) with angular frequency u = eBz/rrie. The operators a) and a 
satisfy the usual bosonic commutation relations and act as raising and lowering operators 
on the system eigenfunctions, respectively. Through the raising operator, we can generate 
the system eigenfunctions in any level n, by starting with the ground state wavefunctions 
n = or the lowest Landau level (LLL). The LLL is characterized by a\l/„=o(-2) = 0, whose 



solutions are given by the normalized wavefunctions 20|, |2l| : 



'^n=0,m{z) = . \ = ( — ) exp ( -ipT Y (3) 

where * denotes complex conjugation, and the quantum number m denotes the angular 
momentum. The degeneracy of the above wavefunction in m implies that one can construct 
general LLL wavefunctions of the form 



= (const)/(z*)exp (^-^^, (4) 

where f{z*) is any arbitrary analytic function of z* . The normalized eigenfunctions for 
arbitrary n and m are given by 

^ = ^— a"^"^n (5) 

^ n,m I — r"- ^0,m- \'-> ) 

vnl 



B. Landau levels with Rashba SOC and Zeeman coupling in the symmetric gauge 

We introduce spin into the system which in the case of 2DEGs in heterostructures enters 
the Hamiltonian through the Zeeman coupling and Rashba SOC {5,6] terms. We assume a 
narrow-gap heterostructure, in which the Rashba SOC term is the dominant contribution, 

n 

while the Dresselhaus SOC term [22)] can be neglected. The Zeeman coupling and the 



Rashba SOC effects are, respectively, described by the matrix operators 1-Lz = gfiBzCz, and 
T^R = (y./h{Jlyax — Ilx<yy), where g is the Lande factor of electrons, fi = eh/ 2m is the Bohr 
magneton, (Ji=x,y,z are the Pauli spin matrices and a is the Rashba SOC parameter. In terms 
of the raising and lowering operators, the Rashba Hamiltonian has the compact form 

«H^^f"/V (6) 

We solve the total Hamiltonian % = Tio + Hr + Tiz for its eigenspinors, ^n,m{z) = 
{^n.m{^)j ^^nmi^))^ J Writing the spinor components as a linear combination of the spin- 
less and normalized eigenfunctions given by Eq. ([5]), 

N / t \ 

n=o \ai J 

where an denotes the up (down) spin coefficient of the nth Landau level, and we use the 
vector notation ^ to denote eigensj»mor solutions. Note that in Eq. (JTj) the summation runs 
over the Landau level index n whilst the angular momentum m is kept constant [s^. The 
Schrodinger equation — i?I)^'7v,m = (I is the 2-by-2 identity matrix) then reads 

(hwia^a + 1) + gf^B, - E) Zto + Zn=o 4^n,m 

-^at + (hwia^a + i) - gfiB, - E) 



0. (8) 



To simplify Eq. (jH]) we utilize the orthogonality of the Landau level wavef unctions, namely 
that n,m\^ n' ,m) = ^n,n' ^oi any valuc of m. Let us denote as M the column vector on the 
left hand side of Eq. ([8]). Now, we consider multiplying both sides of Eq. ([8]) by the state-bra 
{^s,m\, which yields the equation ^* „^(z)M(i "^z = 0, where the integration is performed 
over the entire complex space C. After canceling the orthogonal terms and applying the 
raising and lowering operators, Eq. ([H]) is simplified as 

(nwis + l)+gfiB,~E)al + ^aU,V^\ 
-^al_^^s + {nw{s+\)-g^iB,-E)ai ) 

The resulting equation is a simple system of two equations relating the spinor components 
of state s and its adjacent states s ± 1. Therefore, we can replace s — )■ s — 1 without any loss 
of generality in the top row of Eq. ([9]), to yield a regular eigenvalue equation whose energy 
eigenvalues E are 

E^ = shio± + 2s(a/r)2 (10) 



where ^ = fvjj/2 — g/iBz- In particular, ^ is the energy of the LLL with eigenspinors 
{al_i,al)'^ — (0,1)^. Compared to the energy spectrum of pure Landau levels, the LLL 
energy differs only by the Zeeman term —gfiBz corresponding the electron spins pointing 
antiparallel to the applied magnetic field. Furthermore, in the LLL the wavefunctions do 
not experience any spin splitting from the Zeeman term (since all eigenstates are spin down, 
(7z — —1) and the wavefunctions are completely independent of the Rashba SOC in the 
system. In general, when s 7^ 0, the Zeeman and Rashba SO coupling breaks the spin 
degeneracy and the wavefunctions are highly dependent on the SOC. Let us label the spin- 
split states such that H'^f,^ — Ef'^f^^. The eigenspinor solutions are given by 



= ' , (11) 




where Kg — »«V2s/r ^ j^^± normalization constants. Since the basis wave- 

functions {^s,m} are normahzed, Nf" satisfies \N^\ — 1/ •\/|kJ^"P~+T. Once again, we find 
that the Landau levels are infinitely degenerate since the choice of m does not affect the en- 
ergy eigenvalue. Therefore, taking arbitrary linear combinations in m of the wavefunctions 
yield general solutions as before: 

*? = iV^a„ (12) 

m \ ^ s,m I 

where the normalization constant is determined by the requirement (^^|^^) = 1, i.e. 



C. Gauge invar iance 

We demonstrate gauge invariance of our solutions obtained in the symmetric gauge with 
respect to the wavefunctions in the Landau gauge. The U(l) gauge invariance of electromag- 
netism requires that for a gauge transformation, A ' = A + Vx, the electron wavefunction 
must undergo a corresponding transformation, 

^' = = exp (^-^X^i^, (13) 

in order for the Schrodinger equation to remain invariant in form. In other words the 
electrons acquire an extra phase factor due to the gauge transformation, which implies that 
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physical observables are identical in both gauges. In going from the Landau to the symmetric 
gauge, the required gauge transformation is given by 

Bzxy 

X = ■ (14) 

For simplicity, we illustrate the principle for only the s = 1 eigenstates. Without any loss of 
generality, we can focus on the eigenfunctions that have their cyclotron centres at the origin 
of the system of coordinates, {xq, yo) = 0. Under these set of conditions, the normalized 
eigenfunctions in the Landau gauge have form jl^ . I23I 1 : 

On the other hand, considering Eqs. @, ([5]) and f|TT]) . the general s = 1 wavefunctions in 
the symmetric gauge are of the form 



{2d.,-z/2r') /(^*)exp(-M 



(16) 



Note that Eq. ( |T6|) is obtained after normalizing \l'o(-2) and a^\['o(-z), and substituting the 
explicit expression for the operator. Now, gauge invariance is valid if the same wavefunc- 
tions in the respective gauges are linked via the relation of Eq. ( |T3i) . It therefore suffices 
to construct a wavefunction in the symmetric gauge — via the analytic function f{z*) — for 
which this holds. Consider 

^ fc=0 ^ ^ 

Substituting this choice of / into our symmetric gauge solution, we obtain after some ma- 
nipulation 

where t]'^^'^^ is the normalization coefficient for the up (down) spin branch of the spinor. For 



'^o{z) in the up-spin branch to be correctly normalized, we require rj'^ = a/v^- On the 
other hand, in the down-spin branch we set i]^ = i>/rjj2^/n) to satisfy normalization for 
a^\I'o(-2). This then yields for our symmetric gauge wavefunction 



N I K,f^ \ ( —y'^\ f —ixy 



**(.,y) = ^| |exp(^ exp (19) 



which is just the wavefunction in the Landau gauge multiphed by the gauge transformation 
phase factor: 

= *^(x,y)exp (20) 

III. NUMERICAL SIMULATIONS OF LOCAL SPIN POLARIZATION 

We present some numerical results based on the eigenspinors we derived for the symmetric 
gauge. In Fig.[T]we plot the local, spatial spin polarization {^f^mWi\^t,m) (s = 1; m = 

1; +) level in the symmetric gauge in the 2DEG plane. First of all, we notice that the spin 
distributions are circular in nature, reflecting the spatial probability density distribution of 
the Landau orbits. Of the x, y, 2;-spin components, the most interesting are the in-plane 
X and y components [Figs. [U^a) and (b) respectively] as these components arise from the 
Rashba SOC. For the Rashba Hamiltonian Hr, the effective magnetic field r2(k) is oriented 
in the plane of the 2DEG, and is orthogonal to the in-plane momentum ky, i.e. 0(k) -kn = 0. 
The spin alignment along 0(k) can be seen in Fig. [H if one imagines an electron moving in a 
circular orbit around the origin with a tangential velocity of v = Kk\\/me. Thus, the results 
of our quantum mechanical analysis are in general agreement with the classical picture. 
The 2;-component of the spin, on the other hand, shown in Fig. [T]^c) is uniform along the 
orbit, as it arises from the k- independent Zeeman coupling. For the opposite eigenstate 
(s = l;m = 1;—), the values of the x and ?/-spin components have opposite sign. The 
2;-components, however, are not related by any simple transformation, although the general 
shape of the spatial distributions is the same. Increasing m, one observes an increase in the 
radius of the circular distributions. Finally, the Landau level index s defines the maximum 
number of concentric circular orbits in the electron probability distribution. Therefore, the 
electron states in the higher Landau levels are characterized by a larger number of "ripples" 
in their spin texture. In practice, this spatial modulation of the spin polarization can be 



characterized by means of quantum point contacts (QPC) 2J], as the spatial resolution of 
this technique (~ 100 nm) is comparable to the Landau orbital radii. In particular, one 
could conceive a magnetic-focusing arrangement whereby two QPCs are separated in space 
by twice the cyclotron radius. The first QPC behaves as an electron source, whilst the other 
serves as the collector. Under the influence of the magnetic field, electrons from the source 
follow a semi-circular trajectory in the 2DEG with a cyclotron radius, = hkp/eB {kp 
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is the Fermi wavevector of electrons in the source QPC) and are collected by the detector. 
This technique has been used previously to image the trajectory of cyclotron orbits in 2DEG 



systems in the presence of a vertical magnetic field [25|, |26|. Additionally, we have shown 
that the momentum-dependent SOC field manifests itself as a spatially non-uniform in- 
plane spin-polarization along the electron orbits. This spatial variation may be detected 
experimentall y t hrough the use of magneto-optical techniques such as polarized absorption 



spectroscopy 



27|, magneto-optical Kerr rotation [28| or magneto-reflectivity measurements 



29|. A probe-based QPC technique could also be used to tap the spin-current from the 



system locally (assuming that it does not introduce significant local perturbations to the 
system). By placing the probe at optimal positions corresponding to the peak polarization 
values, we could conceivably draw a highly spin-polarized current, thus implementing an 
efficient spin filtering scheme. 

In summary, we studied the spatial spin polarization texture of Landau levels in the presence 
of Rashba SOC. To do so we solved the wavefunctions for the system in the symmetric 
gauge, demonstrating the gauge invariance of our solutions with previously known solutions 
in the Landau gauge. The two-dimensional analysis of the spatial spin dispersion may 
be important for several reasons (i) the momentum-dependent spin-orbit coupling effect is 
clearly seen in the Landau orbits (see Fig. [1]), and this unifies the quantum mechanical and 
classical pictures, (ii) the theoretically predicted spatial spin distribution may readily be 
verified experimentally using standard magnetic-focusing techniques, and (iii) it may find 
useful applications in spintronics, such as efficient spin filtering devices and sources of spin 
polarized current. 
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FIG. 1: (color online) Local spatial distributip^ of spin polarization of (s = l;m = 1;+) level 
obtained in the symmetric gauge, .^^i\cri\^f^i m=i)- '^^^ spatial coordinates are in arbitrary 



